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ABSTRACT 

We examine the energetics of Coronal Mass Ejections (CMEs) with data from the LASCO 
coronagraphs on SOHO. The LASCO observations provide fairly direct measurements of the mass, 
velocity and dimensions of CMEs. Using these basic measurements, we determine the potential 
and kinetic energies and their evolution for several CMEs that exhibit a flux-rope morphology. 
Assuming flux conservation, we use observations of the magnetic flux in a variety of magnetic 
clouds near the Earth to determine the magnetic flux and magnetic energy in CMEs near the 
Sun. We find that the potential and kinetic energies increase at the expense of the magnetic 
energy as the CME moves out, keeping the total energy roughly constant. This demonstrates 
that flux rope CMEs are magnetically driven. Furthermore, since their total energy is constant, 
the flux rope parts of the CMEs can be considered to be a closed system above ~ 2 Rq. 

Subject headings: Sun: activity — Sun: corona — Sun: magnetic fields — Sun: solar-terrestrial relations 



1. Introduction 

Material ejections are a common phenomenon 
of the solar corona. Since the first observation on 
14 December 1971 (Tousey 1973), several thou- 
sands of CMEs have been seen (Howard et al. 
1985; Kahler 1992; Webb 1992; Hundhausen 1997; 
Gosling 1997). Nevertheless, the mechanisms that 
cause a CME and the forces acting on it during 
its subsequent propagation through the corona are 
largely unknown. Of these two issues, the issue of 
CME propagation through the corona is by far 
more amenable. Past observations have provided 
insufficient coverage of the CME development for 
several reasons: restricted field of view of the coro- 
nagraphs, frequent orbital nights and low sensitiv- 
ity of the instruments. Consequently, past stud- 
ies were largely focused on either the phenomeno- 



logical description and classification of CMEs or 
the measurement of average values for the physi- 
cal properties of the events such as speed, mass, 
kinetic energy (Jackson and Hildner 1978; Howard 
ct al. 1985). The study of the CME energetics, in 
particular, was necessarily restricted to a handful 
of well observed events (Rust et al. 1980; Webb et 
al. 1980). Their analysis revealed the importance 
of the (elusive) magnetic energy and established 
that the potential energy dominates the kinetic 
energy. It was also found that the energy resid- 
ing in shocks, radio continua and other forms of 
radiation was insignificant in comparison to the 
mechanical energy of the ejected material. 

The lessons learned from the past resulted in a 
greatly improved set of instruments; the LASCO 
coronagraphs (Brueckner et al. 1995), aboard the 
SOHO spacecraft (Domingo et al 1995). The loca- 
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tion of the spacecraft at the LI point permits the 
continuous monitoring of the Sun while the combi- 
nation of the three LASCO coronagraphs provides 
an unprecedented field of view from 1.1 R@ to 30 
Rq. The replacement of videcons with CCD de- 
tectors and the very low stray light levels of the 
coronagraphs have led to a vast sensitivity im- 
provement. It is now possible to routinely follow 
the dynamical evolution of a CME. Here, we com- 
pute basic quantities; mass, velocity and geom- 
etry and derive quantities such as the potential, 
kinetic and magnetic energies of CMEs as they 
progress through the outer corona into the helio- 
sphere. To our knowledge, this is the first time 
that detailed observations of the dynamical evolu- 
tion of these quantities has been presented. These 
measurements are expected to provide concrete 
observationally-based constraints on the driving 
forces in CME models. For this study, we focus 
on a group of CMEs that share a common charac- 
teristic; namely, they resemble a helical flux rope 
in the C2 and C3 coronagraph images. We choose 
these events for three reasons: (i) the area of a 
CME that corresponds to the flux rope is usually 
easily identifiable in the coronagraph images, (ii) 
their appearance can be related to the flux rope 
structures measured in-situ from Earth-orbiting 
spacecraft, and (iii) there has been extensive the- 
oretical and observational interest for this class of 
CMEs. 

Several CMEs observed with the LASCO in- 
strument exhibit a helical structure like that of a 
flux rope (Chen ct al. 1997; Dere et al. 1999; Wood 
ct al. 1999). The theoretical basis for flux rope 
configurations in solar and interplanetary plas- 
mas is well established (e.g., Gold 1963; Goldstein 
1983; Chen and Garren 1993; Low 1996; Kumar 
and Rust 1996; Guo et al. 1996; Wu, Guo & Dryer 
1997). These treatments envisage the helical flux 
rope as a magnetic structure that resides in the 
lower corona and erupts to form a CME. There is 
some debate about whether the flux rope is formed 
before the eruption, or whether it is formed as a 
consequence of reconnection processes that lead 
to the eruption. These arguments are related to 
those which consider whether the reconnection oc- 
curs above the sheared arcade which presumably 
forms the flux rope, or below it (Antiochos, Devore 
and Klimchuk 1999). 

Neither the physical mechanisms of the initial 



driving impulse, nor the conditions in the corona 
which determine the subsequent propagation of 
the flux rope are very well known from observa- 
tions. Theoretical models often rely on educated 
guesses to model both the initiation of the CME as 
well as its propagation through the corona. State- 
ments about the energetics, or driving forces be- 
hind CMEs are made on these bases; for instance, 
Chen (1996) and Wu ct al. (1997) show plots of 
the variation of kinetic, potential and magnetic 
energies of CMEs as calculated from their models. 
The measurements we present in this paper are 
expected to yield some clues about the validity of 
the assumptions made in these models. It may be 
emphasized that our measurements are made only 
in the outer corona (2.5 Rq - 30 Rq). They are 
therefore not expected to shed much light on the 
energetics of the flux rope CMEs immediately fol- 
lowing initiation, or on the initiation process itself. 

Our estimates of the magnetic energy of flux- 
rope CMEs are made on the basis of in-situ mea- 
surements of magnetic clouds near the earth. This 
is because flux-rope CMEs ejected from the Sun 
are often expected to evolve into magnetic clouds 
(Rust and Kumar 1994; Kumar and Rust 1996; 
Chen et al. 1997; Gopalswamy et al. 1998). Con- 
versely, in-situ measurements of magnetic clouds 
near the earth suggest that their magnetic field 
configuration resembles a flux rope (Burlaga 1988; 
Lepping, Jones, & Burlaga 1990; Farrugia et al. 
1995; Marubashi 1997). Radio observations of 
moving Type-IV bursts can also probe the mag- 
netic field in CMEs (Stewart 1985; Rust et al. 
1980) but they are so rare that near-Earth mea- 
surements are the most reliable estimates of the 
magnetic flux. It should be borne in mind, how- 
ever, that the precise relationship between CMEs 
and magnetic clouds and the manner in which 
CMEs evolve into magnetic clouds is not very well 
understood (Dryer 1996; Gopalswamy et al. 1998). 
The main reason for this situation is the simple ob- 
servational fact that while CMEs are best observed 
off the solar limb, magnetic clouds are measured 
near the Earth. This issue will hopefully be ad- 
dressed in the near future by the next generation 
of space-borne instruments. 

The rest of the paper is organized as follows: 
We describe our methods of measuring the mass 
and position of a CME and of calculating the dif- 
ferent forms of energy associated with it in § 2. § 3 
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presents the results of our measurements. We dis- 
cuss caveats that accompany these results in § 4 
and draw conclusions in § 5. 

2. Data Analysis 
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Fig. 2. — Solid line: Thomson scattering calcula- 
tion of the angular dependence of the total bright- 
ness of a single electron at a heliocentric distance 
of 10 R Q . The curve is normalized to the bright- 
ness at 0°. Dash-dot line: Ratio of the observed 
relative to the actual mass of a simulated CME, 
centered on the plane of the sky, as a function of 
its angular width (see text). 

2.1. Mass calculations 

White light coronagraphs detect the photo- 
spheric light scattered by the coronal electrons 
and therefore provide a means to measure coro- 
nal density. Transient phenomena, such as CMEs, 
appear as intensity (hence, density) enhancements 
in a sequence of coronagraph images. We com- 
pute the mass for a CME in a manner similar 
to that described by Poland et al. (1981). Af- 
ter the coronagraph images are calibrated in units 
of solar brightness, a suitable pre-event image is 
subtracted from the frames containing the CME. 
The excess number of electrons is simply the ratio 
of the excess observed brightness, B Q i, s , over the 
brightness, B e (6), of a single electron at some an- 
gle, 6, (usually assumed to be 0) from the plane 
of the sky. B e (9) is computed from the Thomson 
scattering function (Billings 1966). The mass, to, 



is then calculated assuming that the ejected mate- 
rial comprises a mix of completely ionized hydro- 
gen and 10% helium. Namely, 

m = iw) ' L9? x 1(r24gr (1) 

After the mass image is obtained, we delineate the 
flux rope by visual inspection, as shown in Fig- 
ure 1. We attempt to circumscribe the cross sec- 
tion of the helical flux rope as seen in the plane of 
the sky. The cavity seen in the white light/mass 
images is taken to be the interior of the flux rope, 
bounded by the helical magnetic field (Figure 1). 
The mass contained in the flux rope is computed 
by summing the masses in the pixels encompassed 
by the flux rope. 

The accuracy of the mass calculations depends 
on three factors: the CME depth and density dis- 
tribution along the line of sight and the angular 
distance of the CME from the plane of the sky. 
All three factors are unknown since the white light 
observations represent only the projection of the 
CME on the plane of the sky. Some additional 
information can be obtained from pB measure- 
ments, but these are only occasionally available. 
Therefore, to convert the observed brightness to a 
mass measurement we have to make an assump- 
tion. Namely, we assume that all the mass in the 
CME is concentrated in the plane of the sky. Since 
CMEs are three-dimensional structures, our calcu- 
lations will tend to underestimate the actual mass. 

To quantify the errors arising from our assump- 
tion, we performed two brightness calculations 
shown in Figure 2. The solid line shows the an- 
gular dependence of the quantity B e {6) in equa- 
tion (1) normalized to its value at 0°. We see that 
our assumption that the ejected mass is always in 
the sky plane (9 = 0°) underestimates the mass 
by about a factor of 2 at angles ~ 50 — 60°. We 
expect that the CMEs in our sample are relatively 
close to the plane of the sky (0 < 50°) since their 
flux rope morphology is clearly visible. 

Next, we investigate the effect of the finite 
width of a CME. We simulate a CME with con- 
stant density per angular bin along the line of 
sight, centered in the plane of the sky at a he- 
liocentric distance of 10 R Q . Using equation (1) 
we calculate the observed mass, m obsi for various 
widths and compare it to the actual mass, m cme 
for the same widths. The dashed line in Figure 2 
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shows the dependence of this ratio, m b s /m cme on 
the width of the CME. For angular widths similar 
to those of the CMEs in our sample (< 60°) the 
mass would be underestimated by about ~ 15%. 

Finally, we estimate the noise in the LASCO 
mass images from histograms of empty sky re- 
gions. The statistics in these areas show a gaussian 
distribution centered at zero, as expected. We de- 
fine the noise level as one standard deviation or 
about 5 x 10 s gr in the C2 frames and 3 x 10 10 gr 
in the C3 frames. The average C2 pixel signal in 
the measured CMEs is 10 times the noise and the 
C2 pixel signal-to-noisc ratio in the mass measure- 
ments is between 10-100. The CMEs get fainter 
as they propagate farther from the sun. There- 
fore, the pixel signal-to-noise ratio in the C3 im- 
ages drops to about 3-4. These figures refer to 
single pixel statistics and demonstrate the quality 
of the LASCO coronagraphs. Our measurements 
are based on statistics of hundreds or thousands 
of pixels for each image. Therefore, the "mass" 
noise in our images is insignificant compared to 
the systematic errors involved in the calculation 
of a CME mass as discussed previously. 

In summary, these calculations suggest that the 
LASCO measurements tend to underestimate the 
CME mass by about 50%, for realistic widths and 
propagation angles. A more detailed analysis of 
CME mass calculations will appear elsewhere. 

2.2. CME Energy calculations 

In this analysis we consider only three forms of 
energy — potential, kinetic, and magnetic energy. 
These energies can be estimated from quantities 
measured directly in the LASCO images like CME 
area, mass and speed. Two of the many other 
forms of energy that can exist in the CME/corona 
system can be estimated based on some assump- 
tions and educated guesses: the CME enthalpy 
U and the thermal energy Et- We will show in 
§ 4 that the thermal energy Et is insignificant. 
There are several uncertainties involved in calcu- 
lating the enthalpy of a CME. Firstly, the tem- 
perature structure of a CME is far from known. 
It is conceivable that is composed of multither- 
mal material. In situ measurements of magnetic 
clouds near the earth reveal a temperature range 
of 10 4 — 10 5 K. Furthermore, it is not clear if the 
gas in the CMEs in the outer corona is in local 
thermodynamic equilibrium. Nonetheless, if we 



assume the CME to be a perfect gas in local ther- 
modynamic equilibrium with equal electron and 
ion temperatures, the enthalpy U can be as large 
as 5Et — bnkT . If we assume a temperature of a 
million degrees K and a mass of 10 15 gr, this yields 
U w 3 x 10 29 ergs. As will be seen later, even this 
upper limit for the enthalpy U is lower than the 
kinetic and potential energies by at least one order 
of magnitude, except in the lower corona where it 
can be comparable to the kinetic energy. Further- 
more, the enthalpy is directly proportional to the 
mass, which, as will be seen later, remains approxi- 
mately constant as the CME propagates outwards. 
We therefore conclude that the enthalpy is a small, 
constant magnitude correction which can be safely 
neglected without affecting the overall conclusions 
regarding CME energetics. 

Potential Energy We define the potential en- 
ergy of the flux rope as the amount of energy re- 
quired to lift its mass from the solar surface. The 
gravitational potential energy is calculated using 
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where and denote the mass and distance 
from sun-center respectively, of each pixel, M Q is 
the mass of the sun, Rq is the solar radius and 
G is the gravitational constant. The summation 
is taken over the pixels comprising the flux rope 
(Figure 1). 

Kinetic Energy We use the center of mass of 
the flux rope to describe its movement. The loca- 
tion of the center of mass relative to the sun center 
is given by 
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where tcm is the radius vector of the center of 
mass and fl is the radius vector for each pixel. The 
summation, as before, is taken over the pixels com- 
prising the flux rope. We calculate tcm for each 
CME frame as it progresses through the LASCO 
field of view. In other words, we compile a table 
of center-of-mass locations versus time, (f C M , t). 
By fitting a second degree polynomial to (feu , t) 
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we obtain the center of mass velocity, vcm and ac- 
celeration acM- The calculation of the speed and 
acceleration as described above has the advantage 
of involving only the measurement of the CME 
center of mass. Once the flux rope is delineated, 
its mass, speed and energetics follow. The kinetic 
energy is simply 

e k = 2 mi "»/ ■ ( 4 ) 

flux rope 

Note that these measurements are based on the 
plane of the sky location of the center of mass. The 
speed used in the calculations is therefore a pro- 
jected quantity and not the true radial speed. It 
follows that the derived kinetic energies are lower 
limits. The same applies for all of our observed 
and derived quantities which facilitates the com- 
parison among the different events. 

Magnetic Energy The calculations of the po- 
tential and kinetic energies of flux rope CMEs are 
made directly from the mass images. On the other 
hand, the values we use for the magnetic energy of 
these CMEs are only estimates because the mag- 
netic field strength in a CME is unknown. In-situ 
measurements by spacecraft like WIND yield the 
magnetic field contained in magnetic clouds ob- 
served near the earth. As mentioned in § 1, helical 
flux-rope CMEs are thought to evolve into mag- 
netic clouds similar to those observed at the earth. 
Therefore, measurements of the magnetic flux con- 
tained in such magnetic clouds are expected to be 
fairly representative of that carried by flux rope 
CMEs. The magnetic energy carried by a flux rope 
CME is defined by 

E M = tt- I B2d V, (5) 

* 7t J flux rope 

where B is the magnetic field carried by the flux 
rope, and the integration is carried out over the 
volume of the flux rope. For a highly conduct- 
ing medium such as the heliosphere, the magnetic 
flux, J BdA, is frozen into the CME as it evolves 
to form a magnetic cloud. The magnetic flux mea- 
sured in-situ is therefore taken to be the same as 
that contained in the CME as it passes through 
the LASCO field of view. We use this frozen flux 
assumption since we feel that it is a simple, phys- 
ically motivated one. Another assumption which 



gives very similar results is conservation of mag- 
netic helicity (Kumar and Rust 1996). The volume 
integral in equation (5) contains another unknown; 
the volume occupied by the flux rope. Assuming a 
cylindrical flux rope with constant magnetic field, 
equation (5) is approximated as 

E M ~±-i(B-Af, (6) 

where A is the area of flux rope as measured in the 
LASCO images and I is the length of flux rope. 
The quantity B ■ A is the magnetic flux frozen into 
the flux rope and is conserved. For our purposes, 
we need, in equation (6), a representative value 
for the magnetic flux of a flux rope. We obtain 
such an estimate from model fits (Lepping, Jones, 
& Burlaga 1990) to several magnetic clouds ob- 
served by WIND between 1995-1998 available at 
http://lepmfi.gsfc.nasa.gov/mfi/mag_cloud_publp.html. 
We only consider clouds that occurred at the same 
time interval as the LASCO CMEs (1997-98). 
From this sample we get the average magnetic flux, 
< B ■ A >= 1.3 ± 1.1 x 10 21 G cm 2 which we put 
in equation (6). The resulting magnetic energy 
uncertainty is then (1.1/1.3) 2 « 70%. To calcu- 
late the magnetic energy, we also need the length 
I of the rope along the line of sight. Since the true 
length of the rope cannot be obtained observa- 
tionally, we assume that the flux rope is expand- 
ing in a self-similar manner, with its length being 
proportional to its heliocentric height; namely, 
I ~ r C M- 

Finally, we emphasize that the magnetic cloud 
data used here are only representative. They are 
not measurements from the same LASCO events 
we analyzed. Also the magnetic flux in individ- 
ual events can differ from the average value we 
adopted. Furthermore, the magnetic field values 
we use refer to the total (toroidal + poloidal) mag- 
netic field contained in the flux rope. The defini- 
tion of B ■ A, however, refers only to the toroidal 
component of the magnetic field which is normal 
to the cross-sectional area of the flux rope. For 
these reasons, it is difficult to ascribe errors to our 
magnetic energy calculations of individual events. 
Therefore, we decided to use the statistical uncer- 
tainty in the average flux to compute the error in 
the magnetic energy which is about 70% as shown 
above. It is unfortunate that the magnetic en- 
ergy measurements are so uncertain and they will 
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Table 1 
CME Event List 



Date 


First C2 appearance 


Position angle 


Angular Width 


Final Speed 




(UT) 


CCW from North (Dcg) 


(Deg) 


(km/sec) 


970223* 


02:55 


89 


63 


920 


970413 c 


16:12 


260 


42 


520 


970430 a 


04:50 


83 


70 


330 


970813 d 


08:26 


272 


36 


350 


971019 b 


04:42 


90 


77 


263 


971030 


18:21 f 


85 


50 


215 


971031 


09:30 


260 


54 


476 


981101 b 


20:11 


272 


57 


264 


980204 


17:02 


284 


43 


420 


980224 


07:55 


88 


32 


490 


980602 e 


09:37 


246 


47 


600 



a Wood et al. (1999) 

b Dere et al. (1999) 

c Chen et al. (1997) 

d Andrcws and Howard (1999) 

°Phmkett, Vourlidas, & Simbcrova (2000) 

f The time refers to the previous day, 97/10/29. 



6 



continue to be so until direct observations of the 
coronal magnetic field become available. 

3. Results 

For our analysis, we searched the LASCO 
database for CMEs with clear flux rope morpholo- 
gies. We picked 11 events for which we compiled 
the evolution of the mass and velocity of the cen- 
ter of mass and the potential, kinetic and mag- 
netic energies as the CME progressed through the 
LASCO C2 and C3 fields of view. For reference 
purposes we present a list of the events in Table 1. 
The information for the 1997 CMEs is taken from 
the LASCO CME list compiled by Chris St. Cyr 
(http: / /lasco-www. nrl.navy.mil/cmelist. html) ex- 
cept for the final speeds in the last column that 
refer to the center of mass of the fluxropes and 
were calculated by us. Further information on 
source regions and associated photospheric/low 
corona emissions for some of these events can be 
found in the references noted in the table. 

Our measurements are shown in Figures 3-6. 
The horizontal axis denotes heliocentric height in 
solar radii. Each row is a separate CME event, la- 
beled by its date of observation by the LASCO/C2 
coronagraph. The left panels show the evolution 
of the potential, kinetic, magnetic and total en- 
ergy in the CME. The total energy is the sum of 
the potential, kinetic and magnetic energies. The 
right panels show the evolution of the flux rope 
mass and the center-of-mass speed. As discussed 
in § 2, a second degree fit to (feu ,t) yields the 
acceleration of the center of mass Scm- The ra- 
dial component of acM is also shown in this panel. 
The dash-dot line, visible in some plots, marks the 
escape speed from the Sun as a function of height. 

An inspection of the plots leads to the follow- 
ing overall conclusions that hold for most of the 
events: 

• The total energy (curves marked with +) is 
relatively constant, to within a factor of 2, 
for the majority of the events despite the 
substantial variation seen in the individual 
energies. This suggests that, for radii be- 
tween approximately 3Rq and 30i?©, the 
flux rope part of these CMEs can be con- 
sidered as an isolated system; i.e., there is 
no additional "driving energy" other than 



the energies we have already taken into ac- 
count (potential and kinetic energies of the 
flux rope, and magnetic energy associated 
with the magnetic field inside the flux rope). 

• We see that the kinetic and, (to a lesser de- 
gree) potential energies increase at the ex- 
pense of the magnetic energy, keeping the 
total energy fairly constant. The decrease 
in magnetic energy is a direct consequence 
of the expansion of the CME. It could imply 
that the untwisting of the flux rope might be 
providing the necessary energy for the out- 
ward propagation of the CME in a steady- 
state situation. 

• The center of mass accelerates for most of 
the events, and the CMEs achieve escape 
velocity at heights of around 8-10 R & , well 
within the LASCO/C3 field of view. 

• The mass in the flux rope remains fairly 
constant for some events (e.g., 97/08/13 or 
97/10/30) while other events (e.g., 97/11/01 
or 98/02/04) exhibit a significant mass in- 
crease in lower heights and tend to a con- 
stant value in the outer corona, above about 
10 — 15 R Q . This observation raises the ques- 
tion: why is pile up of preexisting material 
observed only in some flux rope CMEs? We 
plan to investigate this effect further in the 
future. It would also be interesting to exam- 
ine how the mass increase close to the Sun 
relates to interplanetary "snowplowing" ob- 
servations (Webb, Howard, & Jackson 1996). 

The only notable exception is the event of 
98/06/02 which is also the most massive and its to- 
tal energy increases with distance from the center 
of the sun. This CME is associated with an excep- 
tionally bright prominence which may affect the 
measurements. A detailed analysis of this event 
is presented in Plunkett, Vourlidas, & Simberova 
(2000). 

4. Discussion 

The conclusions of the previous section are 
based on a set of broadband white light corona- 
graph observations. The accuracy of the measure- 
ment of any structure (i.e., CME) in such images 
is inherently restricted by three unknowns: the 
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amount and distribution of the material and the 
extent of the structure along the line of sight. 

We addressed the hrst two problems in § 2 
where we showed that for the case of a uniformly 
filled CME extending ±80 degrees out of the plane 
of the sky, we will measure about 65% of its mass. 
Since the potential and kinetic energies are di- 
rectly proportional to the mass, our measurements 
in Figures 3-6 could be higher by as much as 35%. 
The spatial distribution of the material will also 
affect the visibility of the structures we are trying 
to measure. Because we delineate the area of the 
flux rope by visual inspection, we might not be 
following the same cross section as the structure 
evolves. This might account for some of the vari- 
ability of the energy curves. However, we chose the 
CMEs based on their clear flux rope signatures. 
The measurements involve hundreds or even thou- 
sands of pixels per image and therefore we don't 
expect that the trends seen in the data are affected 
by the slight changes in the visibility of the struc- 
ture. 

The widths along the line of sight of the ob- 
served CMEs are difficult to quantify. There is 
no way to measure this quantity with any instru- 
mentation in existence today. Only the magnetic 
energy depends on the width of the flux rope. In 
§ 2.2, we assumed that the width of the flux rope 
is equal to the height of its center of mass which 
implies that its preeruption length is about a so- 
lar radius. Prominences and loop arcades of this 
length arc not uncommon features on the solar 
surface. 

As described in § 1, flux rope CMEs are ex- 
pected to evolve into magnetic clouds near the 
earth. This is the basis on which we use in-situ 
data to estimate the magnetic energy carried by 
the flux rope CMEs (§ 2). In § 2, we also estimated 
that the overall normalization of the magnetic en- 
ergy curve is uncertain by about 70%. 

In summary, none of the above errors can affect 
the trends of the curves for a given event. Only the 
magnitudes of the various energies could change. 
Finally, some of the variability of the measured 
quantities could be attributed to the intrinsic vari- 
ability of the corona and/or of the CME structure 
itself and cannot be removed without affecting the 
photometry. For this reason, it is rather difficult 
to associate an error estimate to individual mea- 
surements. Therefore, we decided not to include 



any error bars in our figures. 

The analysis of the CME dynamics in Figures 3- 
6 reveals an interesting trend; namely, the total 
energy remains constant. It appears that the flux 
rope part of a CME propagates as a self-contained 
system where the magnetic energy decrease drives 
the dynamical evolution of the system. All the 
necessary energy for the propagation of the CME 
must be injected in the erupting structures dur- 
ing the initial stages of the event. The notions 
that these CMEs are indeed magnetically driven 
and that the thermal energy contribution can be 
ignored arc further reinforced by the magnitude 
of the plasma (3 parameter (Fig. 7). The calcu- 
lations were performed with the assumption that 
the CME material is at a coronal temperature of 
10 6 K. We see that the CMEs have a very small 
(3 (except the events on 98/02/04 and 98/06/02) 
which increases slightly outwards. It appears to 
tend towards a constant value. Such a behavior for 
the plasma (3 parameter was predicted in the flux 
rope model of Kumar and Rust (1996). We also 
find that the potential energy is larger than the 
kinetic energy. These results confirm the conclu- 
sions from earlier Skylab measurements (see Rust 
et al. (1980) for details). 

The relation between the helical structures seen 
in the coronagraph images and eruptive promi- 
nences is still unclear. In our sample, only half 
of the CMEs have clear associations with eruptive 
prominences (e.g., 97/02/23). No helical struc- 
tures are visible in pre-eruption EIT 195A images, 
in agreement with past work (Dere et al. 1999). 
On the other hand, the flux rope of the event on 
98/06/02 is very clearly located above the erupt- 
ing prominence and there is strong evidence that it 
was formed before the eruption (Plunkctt, Vourli- 
das, & Simbcrova 2000). It seems, therefore, likely 
that the process of the formation of the flux rope 
is completed during the early stages of the erup- 
tion at heights below the C2 field of view (< 2 
Rq). Such an investigation, however, is beyond 
the scope of this paper. 

Finally, we turn our attention to the evolution 
of the flux rope shape as a function of height. We 
proceed by comparing the velocity of the CME 
front to its center of mass velocity. Because the 
visual identification of points along the front can 
be influenced by visibility changes as the CME 
evolves, it is susceptible to error. A better method 
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is to use a statistical measure for the location of 
the front such as the center of mass. Hence, the 
location of the front is defined as the center of 
mass of the pixels that lie within 0.1 Rq of the 
front of the flux rope and within ±25° of the radial 
line that connects the sun center with the center 
of mass. The velocity of the front, Vf is calcu- 
lated in the same manner as vqm (§ 2.2). The 
comparison of the two velocity profiles for some 
representative events is shown in Figure 8. Six of 
the eleven CMEs have profiles similar to 97/08/13 
(self-similar expansion) or 97/10/30 (no expan- 
sion) , while five show a progressive flattening such 
as 97/04/13 or 97/11/01, similar to that found in 
Wood et al. (1999). Some theoretical flux rope 
models also predict flattening of the flux rope as 
it propagates outwards (Chen et al. 1997; Wood 
ct al. 1999). 

5. Conclusions 

We have examined, for the first time, the en- 
ergetics of 11 flux rope CMEs as they progress 
through the outer corona into the heliosphere. 
The kinetic and potential energies are computed 
directly from calibrated LASCO C2 and C3 im- 
ages, while the magnetic energy is based on es- 
timates from near-Earth in-situ measurements of 
magnetic clouds. These results are expected to 
provide constraints on flux rope models of CMEs 
and shed light on the mechanisms that drive such 
CMEs. These measurements provide no informa- 
tion about the initial phases of the CME (at radii 
below ~ 2Rq). All the measurements and conclu- 
sions hold for heights in the C2 and C3 fields of 
view; between 3 and 3Oi?0. The salient conclu- 
sions from an examination of 11 CMEs with a flux 
rope morphology are: 

• For relatively slow CMEs, which constitute 
the majority of events, 

— The potential energy is greater than the 
kinetic energy. 

— The magnetic energy advected by the 
flux rope is given up to the potential 
and kinetic energies, keeping the total 
energy roughly constant. In this sense, 
these events are magnetically driven. 

• For the relatively fast CMEs with velocities 
> 600 km/s (97/02/23, 98/06/02), 



— The kinetic energy exceeds the poten- 
tial energy by the time they reach the 
outer corona (above ~ 15i?©). 

— The magnetic energy carried by flux 
rope is significantly below the potential 
and kinetic energies. 



We thank D. Spicer for the initial discussions 
that led to this paper and the referee for his/her 
constructive comments. SOHO is an international 
collaboration between NASA and ESA and is part 
of the International Solar Terrestrial Physics Pro- 
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tium of institutions: the Naval Research Labo- 
ratory (Washington, DC, USA), the University of 
Birmingham (Birmingham, UK), the Max-Planck- 
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Fig. 1. — A typical flux-rope CME seen at the eastern half of the LASCO/C3 coronagraph. The small black 
circle inside the occulter represents the solar disk. The helical structure reminiscent of a flux rope projected 
on the plane of the sky is easily discernible. The dashed circle in the lower left panel demarcates the flux 
rope part of the CME used in our calculations. The location of the center of mass is marked with an 'x' and 
its vector is r CM . 
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Fig. 3. — LASCO measurements of CMEs. Each row is a separate event, labeled by its observation date 
in the C2 coronagraph. The horizontal axis denotes heliocentric height (in R©). Left panels: Evolution of 
potential (dash), kinetic (dash-dot) , magnetic (solid) and total (solid with pluses) CME energies. The total 
energy is the sum of the potential, kinetic and magnclfc energies. Right panels: Evolution of the mass (solid 
line with diamonds) and the center-of-mass speed (dashed line with asterisks). The derived acceleration is 
also shown. The dash-dot line, visible in some plots, is the escape speed from the Sun as a function of height. 
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Fig. 4. — Same as Figure 3. 
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Fig. 5. — Same as Figure 3. 
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Fig. 6. — Same as Figure 3. 
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Fig. 7. — Plot of plasma (3 as a function of height for the flux rope CMEs. We assumed that the CME 
material is at a temperature of 10 6 K. 
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Fig. 8. — Comparison of the speeds between the front (crosses) and the center of mass (asterisks) of the flux 
rope for four representative events. The accelaration resulting for a polynomial fit to these curves is also 
shown where a and aj are the acceleration of the center of mass and the front, respectively. 
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